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ABSTRACT 


This  thesis  forms  part  of  a  larger  study  that  aims  to  develop  a  renewable  energy 
demonstration  plant  at  the  Naval  Postgraduate  School  Turbopropulsion  Laboratory.  The 
architecture  and  design  approach  of  the  demonstration  plant  is  outlined  in  this  thesis. 
While  all  the  components  of  the  system  are  commercially  available,  the  integration  of  the 
components  is  challenging.  The  results  of  the  design  approach  presented  the  optimal  way 
of  integrating  wind  turbines,  an  electrical  system,  chiller  units,  and  thermal  storage  tanks. 
Modular  ice  thennal  tanks  with  polypropylene  tubing  were  chosen  for  storage.  The  ice 
thennal  storage  units  were  selected  over  battery  storage  as  they  are  more  cost  effective 
and  potentially  safer. 

A  statistical  analysis  was  perfonned  using  wind  data  from  Monterey  Airport, 
which  was  beneficial  for  choosing  which  wind  turbines  to  implement  in  the  system.  The 
analysis  determined  that  total  energy  captured  by  two,  4-kW  vertical  axis  wind  turbines 
was  2,554.8  kW-hours  annually.  Additionally,  ANSYS  Fluent  was  used  to  analyze  the 
ice  growth  around  the  tubing  at  various  ice  and  tube  thicknesses.  The  ANSYS  Fluent 
analysis  showed  that  ice  thickness  affects  the  ice  volume  growth  and  change  in  enthalpy 
change  more  than  wall  thickness  affects  these  conditions. 
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I.  INTRODUCTION 


A.  RENEWABLE  ENERGY 

The  primary  sources  of  energy  in  the  United  States  are  currently  being  used  at  an 
unsustainable  rate.  As  a  result,  the  U.S.  government  has  committed  to  increasing  the 
production  and  use  of  renewable  energy  resources.  Diversifying  the  energy  infrastructure 
will  enhance  energy  security  while  greatly  benefiting  the  environment.  Figure  1  shows 
that  36%  of  energy  consumed  in  the  U.S.  in  2011  came  from  petroleum,  with  only  9% 
coming  from  all  renewable  energy  resources  combined. 


Figure  1.  Renewable  Energy  as  Share  of  Total  Primary  Energy  Consumption. 

From  [1]. 
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The  Department  of  Energy  projects  that  by  2040,  23-31%  of  electricity  used  in  the  U.S. 
will  be  from  renewable  resources  [2].  This  marked  increase  indicates  the  importance  of 
developing  renewable  energy  technology  for  the  future. 

1.  Legacy  Renewables 

The  term  “legacy  renewables”  is  intended  to  refer  to  renewable  energy  technology 
that  has  been  extensively  developed  and  used.  This  includes  the  renewable  energy 
sources  that  are  directly  linked  to  the  electrical  grid.  In  the  United  States,  most  large- 
scale  wind  turbine  and  solar  photovoltaic  farms  are  directly  connected  to  the  electrical 
grid.  Both  wind  and  solar  are  intermittent  sources  of  energy.  Wind  power  output  is 
directly  affected  by  wind  conditions  and  air  density,  while  solar  is  affected  by  the  season, 
the  time  of  day,  and  cloud  coverage.  To  maintain  grid  stability,  instantaneous  electrical 
generation  must  match  consumption.  Figure  1  indicates  that  solar  and  wind  power 
combined  only  contributes  1.35%  of  the  total  energy  consumed  in  the  United  States. 
With  these  low  contribution  levels,  maintaining  grid  stability  is  not  difficult.  Wind, 
however,  is  the  fastest  growing  renewable  resource  in  the  United  States  [1].  Grid 
stability  will  be  a  substantial  issue  when  solar  and  wind  resources  become  a  more 
significant  contributor  to  the  electrical  grid.  For  this  reason,  off-grid  solutions,  including 
the  one  proposed  in  this  thesis,  are  becoming  especially  relevant. 

2.  End-Use  Renewables 

“End-use  renewables”  refer  to  renewable  energy  resources  in  self-contained 
systems.  Batteries  are  the  most  common  method  of  storing  off-grid  power;  however, 
thermal  storage  systems  for  heating  and  cooling  applications  are  a  valuable  alternative. 

a.  Batteries 

There  are  a  multitude  of  batteries  used  for  renewable  energy  storage 
systems,  including  nickel-cadmium,  lithium-ion,  alkaline-cell,  and  lead-acid  batteries. 
Lead-acid  batteries  are  the  most  common,  comprising  40-45%  of  batteries  sold 
worldwide  [3].  Their  popularity  is  due  to  their  relatively  low  cost,  coupled  with  a  high 
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electrical  efficiency  and  a  high  power-density.  One  major  disadvantage  of  lead-acid 
batteries  is  the  potential  for  a  hydrogen  evolution  reaction  to  occur,  which  can  cause  the 
battery  to  explode. 

b.  Thermal  Storage 

There  are  numerous  methods  of  thennal  storage;  however,  this  thesis  will 
focus  on  ice-based  thermal  storage.  Heat  of  fusion  refers  to  the  enthalpy  change  that 
occurs  when  a  substance  changes  phase  from  a  solid  to  a  liquid.  Ice  storage  takes 
advantage  of  the  high  heat  of  fusion  of  water  (334.0  kJ/kg). 

There  are  numerous  advantages  to  using  ice  storage.  Compared  to 
batteries,  ice  storage  is  inexpensive.  Table  1  indicates  that  the  purchase  price  of  thennal 
ice  storage  is  approximately  one-quarter  of  that  of  an  alkaline-cell  or  lead-acid  battery 
and  one-tenth  of  that  of  a  lithium-ion  battery. 


Thermal  ice 

Alkaline-cell 

Lead-acid  battery 

Lithium-ion 

storage 

battery 

battery 

$37 

$140 

$150 

$400 

Table  1.  Purchase  Price  of  Various  Storage  Devices  per  kW-hr.  After  [4], 

Thermal  ice  storage  is  significantly  less  hazardous  than  battery  storage.  In  addition,  ice 
storage  does  not  have  a  limiting  life  cycle.  Batteries  can  only  be  recharged  a  given 
number  of  times,  while  the  ice  in  thermal  storage  tanks  can  freeze  and  unfreeze  for  an 
unlimited  number  of  cycles.  Figure  2  shows  a  plot  of  gravimetric  energy  density  versus 
volumetric  energy  density.  This  plot  shows  that  the  weight  to  size  ratio  of  an  ice  storage 
system  is  comparable  to  that  of  a  lead-acid  battery. 
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Figure  2.  Gravimetric  energy  density  versus  volumetric  energy  density.  After  [5], 

B.  MILITARY  APPLICATIONS 

The  Department  of  Defense  is  the  largest  single  consumer  of  energy  in  the  United 
States.  In  FY201 1,  DoD  spent  $17  billion  on  fuel,  approximately  two-thirds  of  which  came 
from  petroleum-based  products  [6].  As  such,  DoD  has  become  proactive  about  establishing  an 
energy-secure  force  by  both  improving  energy  efficiency  and  investing  in  renewable  energy. 
There  are  operational,  economic,  and  enviromnental  benefits  to  doing  so  [7]. 

1.  Operational 

There  are  many  operational  benefits  that  come  with  energy  security.  The  U.S.’s 
dependence  on  foreign  oil  is  inextricably  tied  to  national  security,  and  therefore  military 
strategy  and  procedure.  A  chart  of  conventional  oil  reserves  around  the  globe  is  shown  in 
Figure  3 
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Figure  3.  Proven  conventional  oil  reserves  (2008).  From  [8]. 

Of  these  nations,  only  Mexico  and  Canada  are  considered  “strategically  reliable 
sources  of  energy  to  the  United  States,”  [8].  This  renders  the  U.S.  extremely  vulnerable 
to  the  will  of  strategically  unreliable  foreign  nations.  Achieving  energy  security  would 
lessen  the  dependence  on  these  foreign  nations  while  improving  “the  range,  endurance, 
and  reliability  of  ground,  air,  and  naval  forces,”  [9]. 

One  common  strategic  benefit  of  attaining  energy  security  that  is  frequently  cited 
by  military  officials  is  reducing  the  importance  of  fuel  lines  and  convoys.  A  single 
forward  operating  base  in  Afghanistan  requires  at  least  300  gallons  of  diesel  per  day  [8]. 
Fuel  lines  and  convoys  are  essential  to  providing  fuel  to  these  bases,  which  are  spread 
throughout  the  country.  Convoys  are  extremely  vulnerable  to  IED  attacks,  environmental 
mishaps,  and  traffic  accidents  as  they  travel  from  base  to  base.  In  fact,  there  is  a  direct 
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relationship  between  the  amount  of  fuel  consumed  in  Afghanistan  and  troop  casualties,  as 
shown  in  Figure  4. 

US  total  annual  casualties  (killed  and  wounded)  in  Afganistan 

3.500 
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3,000 

2.500 

2,000 

1.500 

1,431 

1,000 

500 
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0  10  20  30304  40  50  60  656  70 

US  average  monthly  fuel  consumed  in  Afganistan  (Millions  of  gallons) 


FY14 

^ 

\ 

\ 

\ 

\ 

V.  +124%  Increase  ^ 

'  +17.5%  CAGR  ^ 

s 

i 

\ 

\ 

\ 

\ 

FY08 

•  X 

FY09 

FY07 

• 

A- 

^  116%  ) 

FY06 

FY02  X  FY05 
/  *  FY03 

Figure  4.  U.S.  total  annual  casualties  in  OEF  versus  U.S.  fuel  consumed.  From  [8]. 

Reducing  the  dependence  on  petroleum  and  utilizing  other  energy  resources  will 
greatly  simplify  the  logistics  of  an  operation.  With  less  fuel  to  be  moved,  fewer 
casualties  will  occur. 

2.  Economic 

As  discussed  in  the  “Operational”  section,  the  U.S.  is  dependent  on  strategically 
unreliable  foreign  nations  for  fuel  import.  This  makes  the  DoD  vulnerable  to  fluctuations 
in  availability  and  price  of  petroleum.  Diversifying  its  energy  supply  with  renewable 
resources  and  becoming  more  efficient  will  reduce  petroleum  demand  and  increase  DoD 
self-sufficiency  during  this  financially  delicate  time. 

It  is  also  important  to  note  the  “fully-burdened”  cost  of  fuel,  which  accounts  for 
transportation  cost  of  the  fuel,  as  wells  as  the  cost  of  the  air  and  land  forces  used  to 
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protect  it.  The  Marine  Corps  estimated  that  the  “fully-burdened”  cost  of  a  gallon  of  fuel 
ranged  from  $9  to  $16  per  gallon  for  land  deliveries  and  from  $29  to  $31  per  gallon  for 
air  deliveries  [10]. 

3.  Environmental 

While  the  main  goal  of  the  DoD  in  becoming  energy-secure  is  to  save  lives  and 
reduce  spending,  it  is  important  that  military  and  national  policy  are  in  alignment.  The 
DoD’s  current  energy  policy  is  in  agreement  with  the  United  States’  energy  policy,  which 
emphasizes  reducing  greenhouse  gas  emissions  and  reliance  on  fossil  fuels  in  order  to 
protect  the  environment  from  further  damage. 

4.  “Net-zero”  Installations 

The  ultimate  goal  of  the  DoD  and  the  U.S.  Department  of  Energy  is  to  establish 
“net-zero”  energy  installations.  A  “net-zero”  site  is  defined  as  one  that  “produces  as 
much  energy  on-site  from  renewable  energy  generation  or  through  the  on-site  use  of 
renewable  fuels,  as  it  consumes  in  its  buildings,  facilities,  and  fleet  vehicles”  [11].  Wind 
turbines  are  one  source  of  renewable  energy  that  is  being  used  to  meet  this  net-zero  goal. 
Currently,  wind  turbines  are  being  used  to  provide  electricity  at  a  number  of  military 
installations,  including  Guantanamo  Bay  and  San  Clemente  Island  [12].  Typically,  these 
systems  implement  battery  storage,  which,  as  discussed  in  section  I.A.2,  can  be  costly 
and  dangerous.  On  bases  within  the  United  States,  battery  storage  is  a  reasonably  safe 
and  viable  option;  however,  on  forward  operating  bases,  especially  those  within  a 
warzone,  having  exposed  batteries  can  be  extremely  hazardous.  Using  ice  as  a  direct 
thermal  storage  method,  as  this  project  proposes,  has  the  potential  to  save  fuel  and 
therefore  money.  In  addition,  it  may  be  safer  in  forward-operating  areas. 

5.  Energy  System  Technology  Evaluation  Program 

The  Office  of  Naval  Operations  has  provided  funding  to  the  Energy  System 
Technology  Evaluation  Program  (ESTEP).  This  program  supports  both  technology-  and 
policy-based  projects  focusing  on  various  energy-related  topics.  The  projects  use  Navy 
and  Marine  Corps  facilities  and  resources  for  testing.  At  NPS,  numerous  ESTEP-fimded 
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projects  are  conducted  within  in  the  Operations  Research  and  Mechanical  and  Aerospace 
Engineering  Departments.  This  project,  which  focuses  on  the  implementation  of  a 
renewable  energy-based  cooling  system,  was  funded  by  ESTEP. 

C.  PROJECT  OBJECTIVES 

The  ultimate  goal  of  this  project  is  to  develop  a  demonstration  plant  that  uses 
wind  to  power  a  commercial  chiller  unit  at  the  Naval  Postgraduate  School  (NPS) 
Turbopropulsion  Laboratory.  The  excess  cooling  will  be  used  to  freeze  water,  which  will 
be  stored  in  a  thennal  storage  unit.  Wind  power  is  intermittent,  so  when  the  turbines  are 
not  spinning  and  producing  power,  the  ice  storage  unit  will  discharge  to  cool  the  building. 
While  ice  storage  has  become  fairly  common  for  cooling  purposes,  there  has  yet  to  be  a 
system  that  uses  wind  or  another  renewable  energy  source  to  power  such  a  system. 

Ultimately,  this  demonstration  plant  would  serve  as  a  model  for  forward  operating 
bases  and  remote  stations.  The  ice  storage  makes  the  system  grid-independent,  and  in 
addition,  commercial  chillers  can  be  readily  purchased,  transferred  to,  and  set  up  in 
international  areas.  Ideally,  a  system  like  this  could  be  used  to  cool  data  centers  or  other 
spaces  in  forward-operating  bases.  Figure  5  shows  that  space  cooling  represents  13.1% 
of  energy  consumed  in  commercial  spaces  in  2005.  This  indicates  the  significance  of 
developing  renewable  energy  resources  that  caters  to  this  application. 
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Figure  5.  Energy  usage  breakdown  in  commercial  spaces  in  2005.  From  [13]. 


This  system  is  in  direct  accordance  with  the  military’s  goal  to  reduce  petroleum 
usage  and  to  increase  the  use  of  renewable  energy  resources. 

D.  LITERATURE  REVIEW 

Although  research  into  thermal  energy  storage  is  ongoing,  ice  storage  for  cooling 
has  been  implemented  in  various  cases,  as  shown  in  cases  like  [14],  [15],  and  [16].  A 
wind-driven  chilling  system  using  ice  storage  has  yet  to  be  developed.  Holly  Chrystle 
Davis  submitted  a  graduate  thesis  that  explored  wind-electric  ice  making  to  the 
University  of  Colorado  in  1994.  A  variable-frequency  wind  generator  powered  an  ice 
maker,  which  was  used  to  preserve  fish  and  produce  in  developing  nations.  The 
University  of  Colorado,  the  National  Renewable  Energy  Laboratory,  and  Bergey 
Windpower  Company,  Inc.  participated  in  setting  up  this  demonstration  plant.  Davis 
showed  that  a  variable-frequency  wind  turbine  can  operate  a  conventional  ice  maker, 
although  she  experienced  problems  with  “large  start-up  currents  and  corresponding 
voltage  drops  that  occurred”  when  an  ice  maker  was  turned  on  [17].  Other  than  this 
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small-scale  project,  a  renewable  energy  powered  cooling  system  using  ice  thennal 
storage  is  a  novel  idea,  and  as  such,  limited  research  is  available. 
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II.  SYSTEM  ANALYSIS 


The  goal  of  the  system  analysis  was  to  determine  the  total  energy  captured  by  a 
wind  turbine  in  one  year.  To  accomplish  this,  a  statistical  analysis  of  wind  data  was  done 
using  MATLAB. 

A.  STATISTICAL  ANALYSIS  OF  WIND  DATA 

Weather  data  statistics,  including  the  wind  from  the  Monterey  Airport  was 
obtained  from  1980  to  August  2012.  The  data  included  measurements  of  the  wind  speed 
and  direction  at  various  times  throughout  the  day.  Before  1998,  the  wind  data  was 
incomplete,  with  large  gaps  of  missing  data.  Because  the  data  from  1980-1997  was  not 
complete,  it  was  decided  that  the  information  from  January  1998  to  August  2012  would 
be  used  to  complete  the  statistical  analysis.  The  reduced  data  still  yielded  160,706  data 
points,  or  an  average  of  30  points  per  day.  The  wind  rose  in  Figure  6  shows  the 
frequency  of  the  wind  blowing  from  a  given  direction.  At  the  Monterey  Airport,  which  is 
less  than  a  mile  from  the  NPS  Turbopropulsion  Laboratory  as  the  crow  flies,  the 
prevailing  winds  come  from  the  west.  The  average  overall  wind  speed  was  5.77  mph,  or 
2.58  m/s.  The  median  wind  speed  was  6.00  mph,  or  2.68  m/s. 
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Figure  6.  Wind  Rose:  Monterey,  CA  January  1998-August  2012. 

B.  POWER  CURVES 

A  MATLAB  analysis  of  the  wind  data  was  subsequently  completed.  The 
MATLAB  script  can  be  found  in  Appendix  A.  The  first  step  in  the  analysis  was  to  graph 
the  wind  velocity  vs.  the  year  in  MATLAB,  as  shown  in  Figure  7. 
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Figure  7.  Wind  velocity  vs.  year. 

The  wind  velocity  was  plotted  as  a  cumulative  distribution  function  (CDF).  A 
CDF  is  a  statistical  distribution  that  describes  the  probability  that  a  random  variable  is 
less  than  or  equal  to  the  independent  variable  of  the  function.  This  graph,  shown  in 
Figure  8,  serves  as  a  useful  first  step  in  determining  the  wind  characteristics  of  a  given 
site.  Once  the  wind  characteristics  are  determined,  an  appropriate  turbine  can  be  chosen. 
This  process  was  detailed  by  Taylor  [18]. 
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Figure  8.  CDF  plot  of  wind  velocity. 

The  next  step  in  the  statistical  analysis  was  to  determine  what  the  power  output  of 
the  turbine  would  be  in  the  Monterey  Airport  wind  conditions.  The  power  curve  of  the  4 
kW  vertical-axis  wind  turbine,  which  was  provided  by  Urban  Green  Energy,  is  shown  in 
Figure  9. 
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Figure  9.  Power  output  vs.  wind  speed.  From  [19]. 

The  data  from  Figure  9  was  interpolated  with  the  wind  speed  data  from  Monterey  Airport 
in  order  to  determine  the  power  output  of  a  single  wind  turbine  at  the  NPS 
Turbopropulsion  Laboratory.  The  results  are  shown  in  Figure  10,  which  shows  that  peak 
power  is  limited  to  just  over  4  kW. 
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Figure  10.  Power  vs.  time. 

15 


Energy  is  defined  as  the  amount  of  power  used  over  time.  The  direct  relationship 
between  energy  (E)  and  power  (P)  is  shown  in  Equation  1 . 

E  =  \Pdt  (1) 

In  order  to  determine  the  amount  of  energy  captured  annually  from  a  single  turbine, 
cumulative  trapezoidal  integration  was  used.  In  this  method,  the  cumulative  integral  of 
power  was  calculated  with  respect  to  time  using  trapezoidal  integration.  The  total  energy 
captured  by  the  4  kW  turbine  would  have  been  18,748  kW-hours.  Dividing  this  value  by 
the  total  number  of  years  of  data  yields  1,277.4  kW-hours  per  year.  This  value  is  the 
energy  captured  annually  by  one  turbine.  Since  two  turbines  will  be  used  at  the  NPS 
demonstration  plant,  the  total  energy  captured  will  be  2,554.8  kW-hours  annually. 
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III.  PLANT  ARCHITECTURE 


The  ultimate  goal  of  this  project  was  to  develop  a  working  demonstration  plant  of 
a  thermal  storage  unit  in  a  renewable  powered  cooling  system.  A  typical  cooling  system 
uses  a  large  battery  for  energy  storage,  rather  than  thermal  energy.  An  example  of  this 
type  of  schematic  is  shown  in  Figure  1 1 . 


Figure  11.  Typical  cooling  schematic. 
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The  proposed  demonstration  plant  will  operate  according  to  the  schematic  shown 
in  Figure  12. 


Figure  12.  Proposed  demonstration  plant  schematic. 
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The  demonstration  plant  will  be  set-up  at  the  NPS  Turbopropulsion  Laboratory. 
An  overhead  Google  Maps  view  of  the  lab  is  shown  in  Figure  13. 


Figure  13.  Google  Maps  image  of  NPS  Turbopropulsion  Laboratory. 
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The  SolidWorks  schematic  of  the  lab,  with  a  to-scale  version  of  the  electrical 
system,  chiller,  and  thermal  storage  tanks,  is  shown  in  Figure  14.  The  individual 
components  of  the  system  will  be  discussed  below. 
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Figure  14.  SolidWorks  model  of  NPS  Turbopropulsion  Laboratory,  including  electrical 

system,  chiller,  and  thermal  storage  tanks. 


A.  WIND  TURBINE  AND  WINDY  BOY  INVERTERS 

Wind  turbines  use  the  wind’s  energy  to  produce  electricity.  When  wind  moves 
over  a  turbine’s  blades,  it  creates  a  pressure  differential.  This  pressure  differential  causes 
the  blades  to  turn  around  a  rotor,  which  is  attached  to  a  low-speed  shaft.  The  low-speed 
shaft  is  connected  to  a  high-speed  shaft  by  a  gear  box,  which  increases  the  rotational 
speed.  The  high-speed  shaft  drives  the  generator,  which  produces  alternating  current 
(AC)  electricity.  The  turbine  is  connected  to  a  wind-interface  box,  which  rectifies  the 
current,  converting  “wild”  AC  to  direct  current  (DC).  The  wind-interface  box  is  then 
connected  to  a  wind  inverter,  which  converts  the  DC  back  into  useable  50/60  Hz  AC. 
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In  this  demonstration  plant,  two  4-kW  vertical  axis  wind  turbines  will  be  used  to 
power  the  system.  The  turbines,  from  Urban  Green  Energy,  are  shown  in  Figure  15. 


Figure  15.  Urban  Green  Energy  4  kW  Vertical  Axis  Turbine.  From  [19]. 


These  turbines  will  each  be  connected  to  an  SMA  Windy  Boy  inverter,  which  will  yield 
the  usable  AC 

B.  CHILLER  SYSTEM  AND  THERMAL  ENERGY  STORAGE 

In  this  demonstration  plant,  the  chiller  is  powered  by  the  wind  turbines  to  produce 
cooling.  This  system  will  implement  Trane’s  Air-Cooled  Liquid  Chiller  Model  CGAM. 
This  chiller  is  variable-speed,  meaning  that  its  speed  will  be  set  to  match  the  amount  of 
wind  power  entering  the  system.  Large  variable-speed  chillers  are  common,  but  smaller 
units  are  not  yet  widely  available.  This  chiller  has  a  26.4  kW  (7.5  ton)  capacity  and  will 
be  placed  outdoors. 

The  chiller  will  be  connected  to  thermal  ice  storage  tanks.  When  the  wind  turbines 

are  spinning,  the  derived  power  will  be  used  to  operate  the  chiller.  However,  when  wind 

output  is  low,  the  ice  in  the  storage  tank  will  melt.  The  high  energy  that  can  be  absorbed  by 

21 


melting  will  then  be  used  for  cooling  instead  of  for  running  the  chiller.  CALMAC  Thermal 
Energy  Storage  tanks  1045 A,  shown  in  Figure  16,  will  be  used  in  this  system. 


Figure  16.  CAFMAC  Thermal  Energy  Storage  tanks  1045 A.  From  [20]. 


Each  tank  has  a  thermal  storage  capacity  of  144  kWh  (41  ton-hr),  a  volume  of  ice/water 
of  1550  liters  (410  gallons),  and  a  volume  of  coolant  solution  in  the  heat  exchanger  of 
151  liters  (40  gallons)  [20].  The  tank  contains  a  spiral-wound,  polypropylene-tube  heat 
exchanger  surrounded  by  water,  shown  in  Figure  17. 


Figure  17.  Inside  of  a  CAFMAC  Thermal  Energy  Storage  tank  1045 A,  with 

polypropylene  tubing.  From  [20]. 
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C.  STABILITY  BATTERY,  GRID-CONNECTION,  AND  SMA  SUNNY 

ISLAND 

In  a  typical  system,  batteries  are  used  as  the  primary  means  of  storage,  as  shown 
in  Figure  11.  In  this  demonstration  plant,  the  battery  is  required  mainly  for  transient 
stability.  Battery  use  is  only  required  for  30  minutes  at  25%  power  for  lubrication 
immediately  after  startup.  As  previously  discussed,  batteries  are  costly  and  potentially 
dangerous.  They  would  be  eliminated  for  an  ideal  system  on  a  forward-operating  base. 

Capacitors  can  be  used  if  the  battery  system  can  be  designed  small  enough. 
Capacitors  charge  and  discharge  much  faster  than  batteries  and  are  stored  without  charge, 
making  them  ideal  for  transport.  Plus,  a  capacitor  can  be  cycled  many  more  times  than  a 
typical  battery. 

The  demonstration  plant  will  also  be  connected  to  the  grid.  This  is  not  always 
feasible  on  forward-operating  bases,  where  a  generator  would  replace  grid-connection. 

The  chiller  that  will  be  used  in  the  demonstration  plant  is  a  3 -phase  unit.  Each  of 
the  wind  turbine/SMA  Windy  Boy  units  will  be  connected  to  two  of  the  three  phase  legs. 
The  SMA  Sunny  Island  inverters  will  automatically  balance  the  load. 
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IV.  SOLIDIFICATION  MODELING 


A.  INTRODUCTION 

In  a  typical  thermal  storage  tank,  a  coolant  solution  of  25%  ethylene  glycol  and 
75%  water  is  cooled  and  circulated  through  a  heat  exchanger.  In  order  to  simplify  the 
analysis,  pure  ethylene  glycol  was  used  in  the  simulation.  The  coolant,  at  a  temperature 
of  -5.56°C  (22°F),  extracts  heat  from  the  water  until  the  water  freezes.  This  process 
occurs  until  around  95%  of  the  water  has  frozen.  In  order  to  gain  a  more  detailed 
understanding  of  this  process,  an  analysis  of  this  procedure  was  completed  using  ANSYS 
Fluent.  The  goal  of  this  analysis  was  to  determine  how  the  ice  and  tube  thickness 
affected  the  solidification  process.  A  secondary  goal  was  to  use  the  heat  transfer  rate  and 
the  solidification  rate  of  ice  to  determine  the  volume  of  ice  that  formed  around  the  tube. 

The  original  goal  of  this  analysis  was  to  complete  a  transient  analysis  of  the 
solidification  of  ice  around  the  tube.  Unfortunately,  this  analysis  proved  very  difficult  to 
complete,  and  a  steady-state  approach  had  to  be  used. 

B.  MODELING  THE  SYSTEM 

The  process  began  by  modeling  the  tube  in  Design  Modeler.  A  quarter  section  of 
the  tube  was  modeled  in  order  to  simplify  the  simulation  and  minimize  run  times. 
CALMAC  storage  tanks  are  being  used  in  this  design  project.  The  company  provided  the 
dimensions  of  the  tube  [20],  which  have  a  thickness  of  0.0002  m,  diameter  of  0.02m,  and 
length  of  0.2  m.  The  tube  was  then  surrounded  by  a  layer  of  ice  with  a  pre-determined 
thickness.  The  Design  Modeler  sketches  are  shown  in  Figure  18  and  Figure  19. 


25 


001$  0046 


Figure  18.  Quartered  section  of  thermal  ice  storage  system  tubing. 
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Figure  19.  Quartered  section  of  thermal  ice  storage  system  tubing — Detail. 

Once  the  sketch  was  complete  in  Design  Modeler,  it  was  uploaded  into  ANSYS 
Workbench.  Each  component  of  the  sketch,  which  included  the  inner  tube  which 
contained  the  ethylene-glycol  solution,  the  tube  itself,  and  the  outer  ice  layer,  was  meshed 
separately.  Figure  20  shows  the  mesh  of  the  tube,  with  a  thickness  of  0.0002  m  and  an 
ice  layer  with  thickness  of  0.0202  m. 
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Figure  20.  Mesh  of  quartered  tube  section. 

The  Sweep  Method  was  used  in  the  meshing  of  the  ice  layer  in  order  to  account 
for  the  curvature  of  the  tube,  and  the  Inflation  tool  was  used  for  the  mesh  of  the  inner 
tube,  as  shown  clearly  in  Figure  2 1 . 
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Figure  2 1 .  Sweep  and  inflation  methods  of  meshing. 


A  numerical  axisymmetric  analysis  was  subsequently  completed  in  ANSYS 
Fluent.  The  Energy,  Radiation  (Discrete  Ordinates),  and  Viscous  (Realizable  k-epsilon, 
Enhanced  Wall  function)  models  were  all  turned  on  for  the  analysis.  The  inner  tube 
material  was  designated  as  ethylene-glycol  with  an  inner-inlet  temperature  of  -5.65°C 
(21.83°F).  The  tube  was  designated  as  polypropylene.  The  outer  layer  was  designated  as 
ice  with  an  inlet  temperature  of  0.00°C  (32.00°F).  Turbulent  flow  ethylene  glycol  was  set 
to  move  through  the  tube  at  4.4  m/s.  This  speed  was  derived  using  the  formula  for 
Reynolds  number.  Reynolds  number  (Re)  for  pipe  flow  is  defined  according  to  Equation  2. 


Re  = 


p\D 


A 


(2) 


In  Equation  2,  p  is  density,  v  is  average  velocity,  D  is  diameter,  and  p  is  dynamic 
viscosity.  For  turbulent  pipe  flow,  Reynolds  number  must  be  greater  than  4000. 
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Knowing  this  and  the  values  of  p,  DH,  and  p  is  density,  v  is  velocity,  D  is  diameter,  and  p, 
it  was  possible  to  solve  for  velocity. 


kg 


4000  = 


1 1 1 1 .4  — x  v  x  0.02m 
m 


0.0157 


N-s 


m 


Solving  for  v  yields  a  velocity  of  2.83  m/s.  In  order  to  guarantee  that  the  velocity  would 
produce  turbulent  flow,  the  velocity  for  the  simulations  was  increased  to  4.4  m/s.  This 
generated  a  Reynolds  number  of  6230,  well  into  the  turbulent  regime. 

Six  runs  using  this  method  were  completed,  each  using  a  different  ice  thickness 
surrounding  the  tube.  The  entire  process  was  then  repeated  using  a  thinner  wall  tube  with 
a  thickness  of  0.0004  m  in  order  to  detennine  the  effect  of  wall  thickness  on  ice  growth. 
Appendix  B  shows  an  example  of  the  temperature  differentials  throughout  the  tube 
during  an  example  trial. 
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c. 


RESULTS 


The  FLUENT  output  report  is  shown  in  Appendix  C.  After  each  run  was 
completed,  the  mass-weighted  average  enthalpy  (J/kg)  was  recorded  at  the  inlet  and  the 
outlet  of  the  tube.  The  difference  between  these  two  values  was  taken 
and  nondimensionalized  to  account  for  the  length  of  the  tube  and  the  fact  that  the  model 
was  a  quarter  section.  These  values  are  shown  in  Appendix  D,  and  the  results  are  shown 
in  Figure  22. 
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Figure  22.  Nondimensionalized  mass-weighed  enthalpy  vs.  ice  thickness. 

The  results  show  that  as  the  thickness  of  the  ice  increases,  the  enthalpy  difference 
between  the  inlet  and  outlet  of  the  tube  decreased.  The  slight  upturn  in  data  at  an  ice 
thickness  of  0.018  m  could  likely  be  resolved  by  increasing  the  resolution  of  the  mesh. 

The  total  heat  transfer  rate  (W)  was  also  recorded  for  the  inlet  and  outlet  of  the 
tube  after  each  run.  The  difference  between  these  two  values  was  taken  and 
nondimensionalized  in  the  same  manner  as  the  mass-weighted  average  enthalpy.  This 
value  was  divided  by  1  second  in  order  to  determine  an  energy  value,  and  then  divided  by 
the  solidification  rate,  or  the  heat  of  fusion,  of  ice  (334  kJ/kg)  in  order  to  determine  the 
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mass  of  ice  produced.  Finally,  the  volume  of  ice  was  determined  by  dividing  this  value 
by  the  density  of  ice  (916.2  kg/m  ).  The  results  are  shown  in  Figure  23. 


Thickness  of  Ice  (m) 

Figure  23.  Volume  of  ice  produced  vs.  thickness  of  ice. 

The  results  show  that  as  the  thickness  of  the  ice  surrounding  the  tube  increased, 
the  volume  of  the  ice  that  is  produced  decreases.  The  numerical  results  are  shown  in 
Appendix  E. 

Wall  thickness  affects  the  results:  a  thinner  wall  resulted  in  a  greater  enthalpy 
difference  between  the  inlet  and  outlet  therefore  a  greater  volume  of  ice  production.  The 
space  between  tubes  in  the  heat  exchanger  was  estimated  to  be  slightly  greater  than  0.02 
m,  or  slightly  larger  than  the  diameter  of  the  tube.  When  the  thickness  of  the  ice 
approached  the  tube  separation  distance,  ice  growth  was  smaller.  However,  it  should  be 
noted  that  in  this  simulation  only  one  tube  was  modeled,  so  tube  separation  was  not  a 
limiting  factor.  The  main  conclusion  that  could  be  drawn  was  that  the  most  important 
factor  for  ice  enthalpy  and  volume  change  in  this  simulation  was  the  thickness  of  the  ice 
around  the  tube. 
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V.  CONCLUSION 


This  thesis  was  integral  for  the  development  of  the  renewable  energy 
demonstration  plant  at  the  NPS  Turbopropulsion  Laboratory.  Developing  renewable 
energy  technology  is  important  to  both  the  United  States  government  and  military.  The 
work  done  in  this  thesis  is  in  accordance  with  the  goals  of  both  of  these  groups  as  they 
move  toward  energy  independence.  Projects  like  this  will  be  especially  useful  in  the 
transition  to  “net-zero”  installations.  A  scaled-down  version  of  this  system  can  be 
especially  valuable  on  forward-operating  bases,  while  a  scaled-up  version  will  be  useful 
for  cooling  data  centers. 

The  analysis  of  the  weather  statistics  was  useful  in  choosing  an  appropriate 
turbine  and  detennining  the  total  power  output  of  the  turbine.  1,277.4  kW-hours  of 
energy  will  be  captured  by  one  turbine  per  year.  Two  turbines  will  be  used,  meaning 
2,554.8  kW-hours  will  be  captured  annually. 

A  separate  analysis  of  ice  solidification  within  the  thennal  ice  storage  system  was 
also  completed.  These  two  analyses  contributed  to  a  more  in-depth  understanding  of  both 
the  chiller  and  thermal  ice  storage  system  mechanisms.  Ultimately  it  was  found  that  the 
thickness  of  the  ice  had  the  most  significant  impact  on  the  change  in  enthalpy  through  the 
tube  and  the  amount  of  ice  formed  around  the  tube. 

The  various  components  of  the  renewable  energy  plant  and  their  complex 
integration  into  the  system  were  thoroughly  explained.  The  system  successfully 
incorporates  wind  turbines,  an  electrical  system,  chiller  units,  thermal  ice  storage  units,  a 
stability  battery,  and  grid  storage. 
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VI.  RECOMMENDATIONS 


The  next  step  going  forward  is  to  actually  set  up  the  demonstration  plant  and 
ensure  that  the  various  components  can  operate  together.  Once  that  is  completed,  the 
architecture  of  the  demonstration  plant  can  be  improved.  An  extensive  study  of  the 
various  components  will  allow  for  the  optimization  of  the  system.  This  will  predicate  the 
use  of  the  system  in  an  alternative  setting,  such  as  a  forward-operating  base.  Specifically, 
a  study  of  what  size  wind  turbine,  what  size  chiller  unit,  and  how  many  thermal  storage 
tanks  to  use  would  be  advantageous.  In  addition,  making  the  system  truly  grid-  and 
battery-independent  can  be  an  area  of  study. 

A  final  important  step  is  the  implementation  of  solar  panels  into  the  system.  Solar 
panels  are  another  source  of  intennittent,  renewable  energy  that  can  be  used  in  parallel 
with  the  wind  turbines  to  power  the  system.  An  analysis  of  solar  conditions,  similar  to 
the  wind  analysis  completed  in  this  thesis,  would  be  beneficial. 
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APPENDIX  A 


%  Matlab  script  file  to  read  in  wind  data  from  Monterey  Airport 

clear  all 
close  all 

%  Constant 

Data . CutlnSpeed  =  4.0;  %[m/s] 

%  Data  is  read  in  from  file 

[NUM, ~ , ~ ] =xlsread ( ' MRY  hourly_observations  1 998-20 12 . xlsx ') ; 

%  Date  Data  is  separated 
Data. String  =  num2str (NUM ( : , 2 ) )  ; 
string 

Data. Year  =  str2num (Data . String (:, 1 : 4 )) ; 

Data. Month  =  str2num (Data . String ( : , 5 : 6) ) ; 

Data. Day  =  str2num (Data . String ( :  ,  7 : 8  )  )  ; 

Data. Hour  =  str2num (Data . String (:, 9 : 10 )) ; 

%  Windspeed  Data  is  separated 
Data. Wind  =  NUM ( : , 4 ) ; 
mph 

Data. Wind  =  1 . 6*Data . Wind/3 . 6 ; 
m/s 

Data.WindCDF  =  sort (Data .Wind (~isnan (Data .Wind) )) ; 

(nans  excluded) 

%  To  eliminate  NaN 
Data .Wind (isnan (Data .Wind) )  =  0; 

%  CDF  data  (p  is  stats  examples) 

Data.p  =  linspace ( 0 , 1 , size (Data . WindCDF, 1 )) ; 

[x ( : , 1 ) , IA, IC]  =  unique (Data .WindCDF)  ;  %  Keep  only  1  data 

at  each  speed 

y  (  :  , 1)  =  (Data.p (IA) )  ' ; 

[x ( : , 2 ) , IA, IC]  =  unique (Data .WindCDF,  ' first ') ;  %  Keep  only  1  data 

at  each  speed 

y ( : , 2)  =  (Data.p  (IA)  )  '  ; 

%  average  of  the  two 
x  =  mean (x, 2 ) ; 
y  =  mean (y, 2 ) ; 

%  Cut  in  speed  probability  is  found 

Data . CutlnSpeedCDF  =  interpl (x, y, Data . CutlnSpeed) ;  %  Probability 

disp('Cut  in  speed  CDF')  Data . CutlnSpeedCDF 

%  Date  is  converted  to  Matlab  date  numbering  system  Data . DateNumber  = 
datenum (Data . Year , Data . Month, Data . Day, Data . Hour, zeros (size (Data . Hour) ) , 
zeros (size (Data. Hour) ) ) ; 


%  Wind  Velocity  in 
%  Wind  Velocity  in 
%  Not-a-Numbers 


%  Converts  to 

%  Strips  out  years 
%  Strips  out  months 
%  Strips  out  days 
%  Strips  out  hours 
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%  Windspeed  is  plotted  against  dates 
figure (1);  close;  figure (1); 

plot (Data . DateNumber, Data . Wind) ;  grid  on 

datetick('x',10, ' keepticks ' , ' keeplimits ' ) ; 
ylabel ( ' Windvelocity  [m/s]');  xlabel ( ' Year ' ) 

%  Plot  integral  of  windspeed  distribution  figure (2);  close;  figure (2); 
plot (Data . WindCDF, Data . p, ' +-b ') ;  hold  on 
plot (x, y, ' -or ' ) 

plot (Data . Cut InSpeed, Data . CutlnSpeedCDF,  '  ok  '  ) 

plot ( [Data . Cut InSpeed  Data . Cut In Speed] ,  [ 0  Data . CutlnSpeedCDF]  ,  '  -  ok  '  ) 
plot ( [ 0  Data . Cut InSpeed] ,  [Data . CutlnSpeedCDF  Data . CutlnSpeedCDF]  ,  ' -ok ' ) 
ylabel (' Cumalitive  Distribution  Function');  xlabel (' Windvelocity 
[m/s]  '  ) 

%  Data  fitting 

%figure(3);  close;  figure (3); 

%plot (Data . WindCDF, log ( -log ( 1 -Data . p) ) )  ; 

clear  Data. String  NUM  %  Memory  cleanup 

[NUM, ~ , ~ ] =xlsread ( ' Turbine_Specif i cat ions . xlsx ' ) 

WindTurb . WindVel=NUM (  :  ,  1 )  ; 

WindTurb . Power=NUM (  :  ,  2  )  ; 

P=interpl (WindTurb . WindVel , WindTurb . Power,  Data . Wind)  ; 

%plot  of  power  vs.  time 

figure (3) ; close; figure (3) ; 

plot (Data . DateNumber, P) ; grid  on; 

datetick ( ' x ' , 10, ' keepticks ' , ' keeplimits ' ) ; 

ylabel ( ' Power  [ kW] ' ) ;  xlabel ( ' Year ' ) 

%plot  of  energy  vs.  time 

E . total=cumtrapz (Data . DateNumber*24 , P) ; 

figure (4) ; close; figure (4) ; 

plot (Data . DateNumber, E . total ) ; grid  on; 

ylabel (' Energy  [kW-hour]  ' ) ;  xlabel (' Year ') ;  datetick (' x ')  ; 

Data . year s= (Data . DateNumber (end) -Data . DateNumber ( 1 ) ) /365; 

E . annual=E . total (end) /Data. years 
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APPENDIX  B 
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Contours  of  Static  Temperature  (k) 


May  07,2013 

ANSYS  FLUENT  14,0  (3d,  dp,  pbns,  rke) 


Figure  24.  Inlet  with  a  tube  thickness  of  0.0002  m  and  an  ice  thickness  of  0.0204  m 
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Contours  of  Static  Temperature  (k)  May  07, 2013 

ANSYS  FLUENT  14.0  (3d,  dp,  pbns,  rke) 

Figure  25.  Outlet  with  a  tube  thickness  of  0.0002  m  and  an  ice  thickness  of  0.0204  m 
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May  07,  2013 

ANSYS  FLUENT  14.0  (3d,  dp,  pbns,  rke) 


Figure  26.  Side  view  of  a  tube  with  thickness  of  0.0002  m  and  an  ice  thickness  of  0.0204 

m 
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APPENDIX  C 


report 


FLUENT 

Version:  3d,  dp,  pbns,  rke  (3d,  double  precision,  pressure-based,  realizable 
k-epsi Ion) 

Release:  14.0.0 
Title: 

Models 


Model 


Settings 


Space 

Time 

Vi  scous 

wall  Treatment 

Heat  Transfer 

Solidification  and  Melting 

Radiation 

Speci es 

Coupled  Dispersed  Phase 
NOx  Pollutants 
SOx  Pollutants 
Soot 

Mercury  Pollutants 


3D 

Steady 

Realizable  k-epsilon  turbulence  model 
Enhanced  wall  Treatment 
Enabled 
Di  sabled 

Discrete  Ordinate  Model 

Disabled 

Di  sabled 

Di  sabl ed 

Di  sabled 

Disabled 

Oi  sabled 


Material  Properties 


Material:  ice_l  (solid) 


Property 


Density 

cp  (specific  Heat) 

Thermal  Conductivity 
Absorption  Coefficient 
Scattering  coefficient 
Scattering  Phase  Function 
Refractive  Index 


Units 

Method 

Value(s) 

kg/m3 

constant 

916.20001 

j/kg-k 

constant 

2050 

w/m-k 

constant 

2.22 

1/m 

constant 

0 

1/m 

constant 

0 

isotropic 

#f 

constant 

1 

Material:  polypropylene  (solid) 


Property 

Units 

Method 

Value(s) 

Density 

kg/m3 

constant 

946 

Cp  (Specific  Heat) 

j/kg-k 

constant 

1800 

Thermal  conductivity 

w/m-k 

constant 

0.15000001 

Absorption  Coefficient 

1/m 

constant 

0 

Scattering  coefficient 

1/m 

constant 

0 

Scattering  Phase  Function 

isotropic 

#f 

Refractive  Index 

constant 

1 

Material:  ethylene-glycol  (fluid) 


Property 


Units  Method  value (s) 


Density 

Cp  (Specific  Heat) 
Thermal  Conductivity 
Viscosity 
Molecular  Weight 
Absorption  coefficient 


kg/m3 

j/kg-k 

w/m-k 

kg/m-s 

kg/kgmol 

1/m 
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constant 

constant 

constant 

constant 

constant 

constant 


1111.4 

2415 

0.252 

0.0157 

62.0482 

0 
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report 


Scattering  coefficient 

1/m 

constant 

0 

scattering  phase  Function 

isotropic 

#f 

Thermal  Expansion  Coefficient 

l/k 

constant 

0 

Refractive  Index 

constant 

1 

Speed  of  Sound 

m/s 

none 

#f 

Material:  ice  (fluid) 

Property 

Units 

Method 

v« 

Density 

kg/m3 

constant 

915 

Cp  (Specific  Heat) 

18.123364  7.809018  000)  (273.15- 
Thermal  conductivity 

j/kg-k 
273.15:  0  0 
w/m-k 

polynomial 

0  0  0) 
constant 

(200-273.15: 

0.0454 

Viscosity 

kg/m-s 

constant 

1. 72e-05 

Molecular  weight 

kg/kgmol 

constant 

18.01534 

Absorption  coefficient 

1/m 

constant 

0 

Scattering  Coefficient 

1/m 

constant 

0 

Scattering  Phase  Function 

isotropic 

#f 

Thermal  Expansion  coefficient 

l/k 

constant 

0 

Refractive  Index 

constant 

1 

Speed  of  sound 

m/s 

none 

#f 

Material:  aluminum  (solid) 

Property  Units  Method  value(s) 

Density  kq/m3 

Cp  (specific  Heat)  j/kg-k 

Thermal  Conductivity  w/m-k 

Absorption  Coefficient  1/m 

scattering  coefficient  1/m 

Scattering  Phase  Function 
Refractive  index 

cell  zone  conditions 
Zones 


name  id  type 


part-outer_fluid  8  solid 
part-pipe  9  solid 
part-inner_fluid  7  fluid 


Setup  Conditions 
part-outer_fluid 
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constant  2719 

constant  871 

constant  202.4 

constant  0 

constant  0 

isotropic  #f 

constant  1 
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condition 


report 


value 


Material  Name 
Specify  source  terms? 
source  Terms 
specify  fixed  values? 

Fixed  values 

.  #f)  (constant  .  0)  (profile  ))) 

Frame  Motion? 

Relative  To  Cell  Zone 

Reference  Frame  Rotation  speed  (rad/s) 

Reference  Frame  X-Velocity  Of  Zone  (m/s) 

Reference  Frame  Y-Velocity  Of  Zone  (m/s) 

Reference  Frame  Z-Velocity  Of  Zone  (m/s) 

Reference  Frame  X-Origin  of  Rotation-Axis  (m) 

Reference  Frame  Y-origin  of  Rotation-Axis  (m) 

Reference  Frame  Z-Origin  of  Rotation-Axis  (m) 

Reference  Frame  x-Component  of  Rotation-Axis 

Reference  Frame  Y-Component  of  Rotation-Axis 

Reference  Frame  Z-Component  of  Rotation-Axis 

Reference  Frame  user  Defined  zone  Motion  Function 

Mesh  Motion? 

Relative  To  Cell  Zone 

Moving  Mesh  Rotation  Speed  (rad/s) 

Moving  Mesh  x-Velocity  of  zone  (m/s) 

Moving  Mesh  Y-Velocity  Of  Zone  (m/s) 

Moving  Mesh  z-Velocity  of  Zone  (m/s) 

Moving  Mesh  X-Origin  of  Rotation-Axis  (m) 

Moving  Mesh  Y-Origin  of  Rotation-Axis  (m) 

Moving  Mesh  z-Origin  of  Rotation-Axis  (m) 

Moving  Mesh  X-component  of  Rotation-Axis 

Moving  Mesh  Y-Component  of  Rotation-Axis 
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ice_l 

no 

((energy)) 

no 

((temperature  (inactive 
no 
-1 
0 
0 
0 
0 
0 
0 
0 
0 
0 
1 

none 

no 

-1 

0 
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report 

Moving  Mesh  Z-Component  of  Rotation-Axis  1 

Moving  Mesh  user  Defined  Zone  Motion  Function  none 

Participates  in  radiation  yes 

Deactivated  Thread  no 

part -pipe 

Condition  Value 


Material  Name 
Specify  source  terms? 

Source  Terms 
Specify  fixed  values? 

Fixed  values 

.  #f)  (constant  .  0)  (profile  ))) 

Frame  Motion? 

Relative  To  Cell  Zone 

Reference  Frame  Rotation  Speed  (rad/s) 

Reference  Frame  x-velocity  of  zone  (m/s) 

Reference  Frame  Y-velocity  Of  zone  (m/s) 

Reference  Frame  Z-Velocity  of  zone  (m/s) 

Reference  Frame  X-origin  of  Rotation-Axis  (m) 

Reference  Frame  Y-origin  of  Rotation-Axis  (m) 

Reference  Frame  z-Origin  of  Rotation-Axis  (m) 

Reference  Frame  X-component  of  Rotation-Axis 

Reference  Frame  Y-Component  of  Rotation-Axis 

Reference  Frame  Z-Component  of  Rotation-Axis 

Reference  Frame  User  Defined  Zone  Motion  Function 

Mesh  Motion? 

Relative  To  cell  zone 

Moving  Mesh  Rotation  Speed  (rad/s) 

Moving  Mesh  X-Velocity  Of  Zone  (m/s) 

Moving  Mesh  Y-Velocity  Of  Zone  (m/s) 

Moving  Mesh  Z-Velocity  Of  zone  (m/s) 
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polypropylene 

no 

((energy)) 

no 

((temperature  (inactive 
no 
-1 
0 
0 
0 
0 
0 
0 
0 

0 

0 

1 

none 

no 

-1 

0 

0 

0 

0 
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report 

Moving  Mesh  X-Origin  of  Rotation-Axis  (m)  0 

Moving  Mesh  Y-Origin  of  Rotation-Axis  (m)  0 

Moving  Mesh  z-Origin  of  Rotation-Axis  (in)  0 

Moving  Mesh  X-Component  of  Rotation-Axis  0 

Moving  Mesh  Y-Component  of  Rotation-Axis  0 

Moving  Mesh  z-Component  of  Rotation-Axis  1 

Moving  Mesh  user  Defined  Zone  Motion  Function  none 

Participates  in  radiation  no 

Deactivated  Thread  no 

part-inner_fluid 

Condition  Value 


Material  Name 


ethylene-glycol 


specify  source  terms? 


no 


Source  Terms  ((mass) 

(x-momentum)  (y-momentum)  (z-momentum)  (k)  (epsilon)  (energy)) 


Specify  fixpri  values7 


no 


Local  Coordinate  System  for  Fixed  Velocities 


no 


Fixed  values  ((x-velocity 

(inactive  .  #f)  (constant  .  0)  (profile  ))  (y-velocity  (inactive  .  #f)  (constant  . 
0)  (profile  ))  (z-velocity  (inactive  .  #f)  (constant  .  0)  (profile  ))  (k  (inactive 
.  # f)  (constant  .  0)  (profile  ))  (epsilon  (inactive  .  #f)  (constant  .  0)  (profile 
))  (temperature  (inactive  .  #f)  (constant  .  0)  (profile  ))) 
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Frame  Motion?  no 

Relative  to  cell  zone  -1 

Reference  Frame  Rotation  speed  (rad/s)  0 

Reference  Frame  X-Velocity  of  Zone  (m/s)  0 

Reference  Frame  Y-Velocity  of  Zone  (m/s)  0 

Reference  Frame  Z-Velocity  of  zone  (m/s)  0 

Reference  Frame  X-Origin  of  Rotation-Axis  (m)  0 

Reference  Frame  Y-origin  of  Rotation-Axis  (m)  0 

Reference  Frame  Z-Origin  of  Rotation-Axis  (m)  0 

Reference  Frame  X-Component  of  Rotation-Axis  0 

Reference  Frame  Y-Component  of  Rotation-Axis  0 

Reference  Frame  z-component  of  Rotation-Axis  1 

Reference  Frame  user  Defined  zone  Motion  Function  none 
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report 


Mesh  Motion?  no 

Relative  To  cell  zone  -1 

Moving  Mesh  Rotation  Speed  (rad/s)  0 

Movinq  Mesh  x-velocity  of  Zone  (m/s)  0 

Moving  Mesh  Y-velocity  of  zone  (m/s)  0 

Moving  Mesh  z-velocity  of  zone  (m/s)  0 

Moving  Mesh  x-Origin  of  Rotation-Axis  (m)  0 

Moving  Mesh  Y-Origin  of  Rotation-Axis  (m)  0 

Moving  Mesh  Z-Origin  of  Rotation-Axis  (m)  0 

Moving  Mesh  X-Component  of  Rotation-Axi s  0 

Moving  Mesh  Y-component  of  Rotation-Axis  0 

Moving  Mesh  Z-Component  of  Rotation-Axis  1 

Moving  Mesh  User  Defined  Zone  Motion  Function  none 
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participates  in  radiation  yes 

Deactivated  Thread  no 

Laminar  zone?  no 

Set  Turbulent  viscosity  to  zero  within  laminar  zone?  yes 

Embedded  subgrid-scale  Model  0 

Momentum  Spatial  Discretization  0 

Cwale  0.325 

cs  0.1 

Porous  zone?  no 

Conical  porous  zone?  no 

X-Component  of  Direction-1  Vector  1 

Y-Component  of  Direction-1  vector  0 
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Z-Cuiriponenl  of  Direction- 1  Vector  0 

X  Component  of  Direction  2  vector  0 

V-Component  of  Direction-2  Vector  1 

Z-Component  of  Direction-2  Vector  0 

X-Component  of  Cone  Axis  vector  1 

Y-Component  of  cone  Axis  vector  0 

Z-Component  of  cone  Axis  vector  0 

x-Coordinate  of  Point  on  cone  Axis  (m)  1 

Y-coordinate  of  Point  on  cone  Axis  0 

z-coordinate  of  Point  on  cone  Axis  (m)  0 

Half  Angle  of  Cone  Relative  to  its  Axis  Cdeg>  0 

Relative  Velocity  Resistance  Formulation?  yes 

Direction-1  Viscous  Resistance  (l/m2)  0 
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report 


Direction-2  Viscous  Resistance  (l/m2) 

Direction-3  viscous  Resistance  (l/m2) 

Choose  alternative  formulation  for  inprrial  resistance? 

Direction-1  Inertial  Resistance  (1/irO 

Direction-2  inertial  Resistance  ( 1/m ) 

Direction-3  inertial  Resistance  (l/m) 

CO  Coefficient  for  Power-Law 

Cl  Coefficient  for  Power-Law 

Porosity 

Equilibrium  Thermal  Model  (if  no,  Non-Equilibrium)? 

Non-Equilibrium  Thermal  Model? 

Solid  Material  Name 
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0 

0 

no 

0 

0 

0 

0 

0 

1 

yes 

no 

aluminum 
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Interfacial  Area  Density  (1/m) 


1 


Heat  Transfer  coefficient  (w/m2-k) 


1 


Boundary  Conditions 


Zones 


name 

id 

type 

outer_outlet 

15 

wall 

outer  inlet 

13 

wall 

wall -part-inner_fluid-part-pipe-shadow 

23 

wall 

wall -part-pipe 

1 

wall 

wal 1 -part-i nner_fl ui d-part-pi  pe 

3 

wall 

inner_inlet 

12 

velocity-inlet 

inner_outlet 

14 

pressure-outlet 

symmetry_l-part-inner_fluid 

17 

symmetry 

symmetry_l-part-pipe 

19 

symmetry 

symmetry_l-part-outer_fluid 

18 

symmetry 

symmetry_2-part-inner_fluid 

20 

symmetry 

symmetry_2-part-pipe 

22 

symmetry 

symmetry_2-part-outer_fluid 

21 

symmetry 

wall 

10046 

wall 

Setup  Conditions 
outer_outlet 
Condition 


value 


wall  Thickness  (m)  0 

Heat  Generation  Rate  (w/m3)  0 

Material  Name  ice_l 

Thermal  bc  Type  1 

Temperature  (k)  300 

Heat  Flux  (w/m2)  0 

Convective  Heat  Transfer  Coefficient  (w/m2-k)  0 

Free  Stream  Temperature  (k)  300 

Enable  shell  conduction?  no 

Wall  Motion  0 

Shear  Boundary  condition  0 
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Define  wall  motion  relative  to  adjacent  cell  zone? 

Apply  a  rotational  velocity  to  this  wall? 

Velocity  Magnitude  (m/s) 

X-Component  of  Wall  Translation 
Y-Component  of  wall  Translation 
z-Component  of  wall  Translation 
Define  wall  velocity  components? 

X-Component  of  wall  Translation  (m/s) 

Y-component  of  wall  Translation  (m/s) 

Z-component  of  wall  Translation  (m/s) 

Internal  Emissivity 
External  Emissivity 
External  Radiation  Temperature  (k) 

Radiation  BC  Type 

X-Component  of  Radiation  Direction 

Y-Component  of  Radiation  Direction 

Z-Component  of  Radiation  Direction 

Theta  Width  of  Beam  (deg) 

9.9923854e-07 

Phi  Width  of  Beam  (deg) 

9.9923854e-07 

(((constant  .  0)  (profile  ))) 

(((constant  .  0)  (profile  ))) 

Apply  Direct  Irradiation  Parallel  to  the  Beam? 

Use  Beam  Direction  from  Solar  Load  Model  Settings 

use  Direct  and  Diffuse  Irradiation  from  Solar  Load  Model  Settings 

Rotation  speed  (rad/s) 

X- Position  of  Rotation -Axis  Origin  (m) 

Y-Position  of  Rotation-Axis  origin  (m) 

Z-Position  of  Rotation-Axis  Origin  (m) 

X-component  of  Rotation-Axis  Direction 
Y-Component  of  Rotation-Axis  Direction 
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yes 

no 

0 

1 

0 

0 

no 

0 

0 

0 

1 

1 

300 

3 

1 

0 

0 


(1) 

no 

no 

no 

0 

0 

0 

0 

1 

0 
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z  component  of  Rotation  Axis  Direction  0 

X-component  of  shear  stress  (pascal)  0 

Y-component  of  shear  stress  (pascal)  0 

Z-component  of  shear  stress  (pascal)  0 

Fslip  constant  0 

Eslip  constant  0 

Surface  tension  gradient  (n/m-k)  0 

Specularity  Coefficient  0 

outer_inlet 

condition  value 


wall  Thickness  (m)  0 

Heat  Generation  Rate  (w/m3)  0 

Mate  ri al  Name  i ce_l 

Thermal  BC  Type  1 

Temperature  (k)  273 

Heat  Flux  (w/m2)  0 

convective  Heat  Transfer  coefficient  (w/m2-k)  0 

Free  Stream  Temperature  (k)  300 

Enable  shell  conduction?  no 

wall  Motion  0 

Shear  Boundary  condition  0 

Define  wall  motion  relative  to  adjacent  cell  zone?  yes 

Apply  a  rotational  velocity  to  this  wall?  no 

Velocity  Magnitude  (m/s)  0.001 

x-component  of  wall  Translation  1 

Y-component  of  wall  Translation  0 

z-component  of  wall  Translation  0 

Define  wall  velocity  components?  no 

X-Component  of  Wall  Translation  (m/s)  0 

Y-Component  of  wall  Translation  (m/s)  0 
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Z-component  of  wall  Translation  (m/s)  0 

Internal  Emissivity  1 

External  Emissivity  1 

External  Radiation  Temperature  (k)  300 

Radiation  BC  Type  3 

x-Component  of  Radiation  Direction  1 

Y-component  of  Radiation  Direction  0 

Z-Component  of  Radiation  Direction  0 


Theta  Width  of  Beam  (deg) 
9. 9923854e-07 

Phi  Width  of  Beam  (deg) 
9.99238S4e-07 

(((constant  .  0)  (profile  ))) 
(((constant  .  0)  (profile  ))) 


Apply  Direct  Irradiation  Parallel  to  the  Beam?  no 

Use  Beam  Direction  from  Solar  Load  Model  Settings  no 

Use  Direct  and  Diffuse  irradiation  from  Solar  Load  Model  Settings  no 
Rotation  Speed  (rad/s)  0 

X-Position  of  Rotation-Axis  Origin  (m)  0 

Y-Position  of  Rotation-Axis  Origin  (m)  0 

Z-Position  of  Rotation-Axis  Origin  (m)  0 

X-Component  of  Rotation-Axis  Direction  1 

Y-Component  of  Rotation-Axis  Direction  0 

Z-Component  of  Rotation-Axis  Direction  0 

X-component  of  shear  stress  (pascal)  0 

Y-component  of  shear  stress  (pascal)  0 

Z-component  of  shear  stress  (pascal)  0 

Fslip  constant  0 

Eslip  constant  0 

Surface  tension  gradient  (n/m-k)  0 

Specularity  Coefficient  0 

wall-part-inner_fluid-part-pipe- shadow 
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Condition  value 


wall  Thickness  (m)  0 

Heat  Generation  Rate  (w/m3)  0 

Material  Name 

al umi num 

Thermal  BC  Type  3 

Temperature  (k)  300 

Heat  Flux  (w/m2)  0 

Convective  Heat  Transfer  Coefficient  (w/m2-k)  0 

Free  Stream  Temperature  (k)  300 

Enable  shell  conduction?  no 

Wall  Motion  0 

Shear  Boundary  Condition  0 

Define  wall  motion  relative  to  adjacent  cell  zone?  yes 

Apply  a  rotational  velocity  to  this  wall?  no 

Velocity  Magnitude  (m/s)  0 

X-Component  of  Wall  Translation  1 

Y-Component  of  Wall  Translation  0 

Z-Component  of  Wall  Translation  0 

Define  wall  velocity  components?  no 

X-Component  of  Wall  Translation  (m/s)  0 

Y-component  of  Wall  Translation  (m/s)  0 

Z-Component  of  Wall  Translation  (m/s)  0 

internal  Emissivity  1 

External  Emissivity  1 

External  Radiation  Temperature  (k)  300 

Radiation  BC  Type  3 

x-component  of  Radiation  Direction  1 

Y-Component  of  Radiation  Direction  0 

z-component  of  Radiation  Direction  0 

Theta  Width  of  Beam  (deg) 
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9.9923854e-07 

Phi  width  of  Beam  (deg) 

9.9923854e-07 

(((constant  .  0)  (profile  ))) 

(((constant  .  0)  (profile  ))) 

(1) 

Apply  Direct  Irradiation  Parallel  to  the  Beam?  yes 

use  Beam  Direction  from  solar  Load  Model  Settings  no 

Use  Direct  and  Diffuse  Irradiation  from  Solar  Load  Model  Settings  no 
Rotation  Speed  (rad/s)  0 

X-Position  of  Rotation-Axis  Origin  (m)  0 

Y-Position  of  Rotation-Axis  Origin  (m)  0 

Z-Position  of  Rotation-Axis  Origin  (m)  0 

X-Component  of  Rotation -Axis  Direction  0 

Y-component  of  Rotation-Axis  Direction  0 

z-component  of  Rotation-Axis  Direction  1 

X-component  of  shear  stress  (pascal)  0 

Y-component  of  shear  stress  (pascal)  0 

Z-component  of  shear  stress  (pascal)  0 

Fslip  constant  0 

Eslip  constant  0 

Surface  tension  gradient  (n/m-k)  0 

Specularity  Coefficient  0 

wall -part-pipe 

Condition  Value 


wall  Thickness  (m)  0 

Heat  Generation  Rate  (w/m3)  0 

Material  Name 

alumi  num 

Thermal  BC  Type  1 

Temperature  (k)  300 

Heat  Flux  (w/m2)  0 
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Convective  Heat  Transfer  coefficient  (w/m2-k)  0 

Free  Stream  Temperature  (k)  300 

Enable  shell  conduction?  no 

Wall  Motion  0 

Shear  Boundary  Condition  0 

Define  wall  motion  relative  to  adjacent  cell  zone?  yes 

Apply  a  rotational  velocity  to  this  wall?  no 

velocity  Magnitude  (m/s)  0 

x-component  of  wall  Translation  1 

Y-component  of  Wall  Translation  0 

Z-Component  of  wall  Translation  0 

Define  wall  velocity  components?  no 

X-Component  of  Wall  Translation  (m/s)  0 

Y-Component  of  Wall  Translation  (m/s)  0 

Z-Component  of  Wall  Translation  (m/s)  0 

Internal  Emissivity  1 

External  Emissivity  1 

External  Radiation  Temperature  (k)  300 

Radiation  BC  Type  3 

X-Component  of  Radiation  Direction  1 

Y-component  of  Radiation  Direction  0 

Z-Component  of  Radiation  Direction  0 

Theta  width  of  Beam  (deg) 

9.9923854e-07 

Phi  Width  of  Beam  (deg) 

9.9923854e-07 

(((constant  .  0)  (profile  ))) 

(((constant  .  0)  (profile  ))) 

(1) 

Apply  Direct  irradiation  parallel  to  the  Beam?  yes 

Use  Beam  Direction  from  solar  toad  Model  settings  no 

use  Direct  and  Diffuse  Irradiation  from  Solar  Load  Model  Settings  no 

Rotation  Speed  (rad/s)  0 

X-Position  of  Rotation-Axis  Origin  (m)  0 
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Y-Position  of  Rotation-Axis  Origin  (m)  0 

Z-Position  of  Rotation-Axis  origin  (m)  0 

x-coniponent  of  Rotation-Axis  Direction  0 

Y-component  of  Rotation-Axis  Direction  0 

Z-Component  of  Rotation-Axis  Direction  1 

X-component  of  shear  stress  (pascal)  0 

Y-component  of  shear  stress  (pascal)  0 

Z-component  of  shear  stress  (pascal)  0 

Fslip  constant  0 

Eslip  constant  0 

Surface  tension  gradient  (n/m-k)  0 

Specularity  coefficient  0 

wall -pa rt-inner_fluid-part-pipe 

Condition  Value 


Wall  Thickness  (m)  0 

Heat  Generation  Rate  (w/m3)  0 

Material  Name 

aluminum 

Thermal  BC  Type  3 

Temperature  (k)  300 

Heat  Flux  (w/m2)  0 

Convective  Heat  Transfer  coefficient  (w/m2-k)  0 

Free  stream  Temperature  (k)  300 

Enable  shell  conduction?  no 

Wall  Motion  0 

Shear  Boundary  Condition  0 

Define  wall  motion  relative  to  adjacent  cell  zone?  yes 

Apply  a  rotational  velocity  to  this  wall?  no 

Velocity  Magnitude  (m/s)  0 

x-component  of  wall  Translation  1 
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Y-component  of  wall  Translation  0 

Z-component  of  wall  Translation  0 

Define  wall  velocity  components?  no 

X-Component  of  Wall  Translation  (m/s)  0 

Y-Component  of  Wall  Translation  (m/s)  0 

Z-Component  of  Wall  Translation  (m/s)  0 

Internal  Emissivity  1 

External  Emissivity  1 

External  Radiation  Temperature  (k)  300 

Radiation  bc  Type  3 

X-Component  of  Radiation  Direction  1 

Y-Component  of  Radiation  Direction  0 

Z-Component  of  Radiation  Direction  0 

Theta  width  of  Beam  (deg) 

9.9923854e-07 

Phi  Width  of  Beam  (deg) 

9.9923854e-07 


(((constant  .  0)  (profile  ))) 

(((constant  .  0)  (profile  ))) 

(1) 

Apply  Direct  Irradiation  Parallel  to  the  Beam?  yes 

Use  Beam  Direction  from  Solar  Load  Model  Settings  no 

Use  Direct  and  Diffuse  Irradiation  from  Solar  Load  Model  Settings  no 

Rotation  Speed  (rad/s)  0 

x-position  of  Rotation-Axis  origin  (m)  0 

Y-Position  of  Rotation-Axis  Origin  (m)  0 

z-position  of  Rotation-Axis  origin  (m)  0 

X-Component  of  Rotation-Axis  Direction  0 

Y-Component  of  Rotation-Axis  Direction  0 

Z-Component  of  Rotation-Axis  Direction  1 

x-component  of  shear  stress  (pascal)  0 

Y-component  of  shear  stress  (pascal)  0 

Z-component  of  shear  stress  (pascal)  0 

Fslip  constant  0 
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Eslip  constant  0 
Surface  tension  gradient  (n/m-k)  U 
Specularity  Coefficient  0 


inner_iniet 


Condition  Value 


Velocity  Specification  Method  2 

Reference  Frame  0 

velocity  Magnitude  (m/s)  4.4 

Supersonic/Initial  Gauge  Pressure  (pascal)  0 
Coordinate  System  0 

x-velocity  (m/s)  0 

Y-velocity  (m/s)  0 

Z-Velocity  (m/s)  0 

X-Component  of  Flow  Direction  1 

Y-Component  of  Flow  Direction  0 

Z-Component  of  Flow  Direction  0 

X-Component  of  Axis  Direction  1 

Y-Component  of  Axis  Direction  0 

Z-Component  of  Axis  Direction  0 

x-coordinate  of  Axis  Origin  (m)  0 

Y-Coordinate  of  Axis  Origin  (m)  0 

Z-Coordinate  of  Axis  Origin  (m)  0 

Angular  velocity  (rad/s)  0 

Temperature  (k)  267.5 

Turbulent  Specification  Method  0 

Turbulent  Kinetic  Energy  (m2/s2)  0.01 

Turbulent  Dissipation  Rate  (m2/s3)  0.1 

Turbulent  intensity  (%)  10 

Turbulent  Length  Scale  (m)  1 

Hydraulic  Diameter  (m)  1 

Turbulent  Viscosity  Ratio  10 

External  Black  Body  Temperature  Method  0 

Black  Body  Temperature  (k)  300 

internal  Emissivity  1 

is  zone  used  in  mixing-plane  model?  no 

inner_outlet 

Condition  Value 


Gauge  Pressure  (pascal)  0 

Backflow  Total  Temperature  (k)  300 

Backflow  Direction  Specification  Method  1 

coordinate  system  0 

X-Component  of  Flow  Direction  1 

Y-component  of  Flow  Direction  0 

Z-Component  of  Flow  Direction  0 

X-Component  of  Axis  Direction  1 

Y-component  of  Axis  Direction  0 

Z-Component  of  Axis  Direction  0 

X-Coordinate  of  Axis  Origin  (ml  0 

Y-Coordinate  of  Axis  Origin  (m)  0 

Z-Coordinate  of  Axis  Origin  (m)  0 

Turbulent  specification  Method  0 

Backflow  Turbulent  Kinetic  Energy  (m2/s2)  1 

Backflow  Turbulent  Dissipation  Rate  (m2/s3)  1 
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Backflow  Turbulent  Intensity  (%)  10 

Backflow  Turbulent  Length  Scale  (m)  1 

Backflow  Hydraulic  Diameter  (m)  1 

Backflow  Turbulent  Viscosity  Ratio  10 

External  Black  Body  Temperature  Method  0 

Black  Body  Temperature  (k)  300 

internal  Emissivity  1 

is  zone  used  in  mixing-plane  model?  no 

Radial  Equilibrium  Pressure  Distribution  no 

Specify  Average  Pressure  Specification  no 

Specify  targeted  mass  flow  rate  no 

Targeted  mass  flow  (kg/s)  1 

Upper  Limit  of  Absolute  Pressure  value  (pascal)  5000000 

Lower  Limit  of  Absolute  Pressure  value  (pascal)  1 


symmet  ry_l-part-i  nner_fl  ui  d 
Condition  Value 


symmetry_l-part-pi pe 
Condition  Value 

symmet  ry_l-part-outer_f 1 ui  d 
Condition  Value 


symmet ry_2-part-i nner_fl ui d 
Condition  Value 


symmet ry_2- part -pipe 
Condition  value 

symmet ry_2-part-outer_fluid 
Condition  value 

wall 

Condition  Value 


wall  Thickness  (m) 

0 

Heat  Generation  Rate  (w/m3) 

0 

Material  Name 

ice_l 

Thermal  BC  Type 

0 

Temperature  (k) 

273 

Heat  Flux  (w/m2) 

0 
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Convective  Heat  Transfer  Coefficient  (w/m2-k)  0 

Free  Stream  Temperature  (k)  300 

Enable  shell  conduction?  no 

Mall  Motion  0 

Shear  Boundary  Condition  0 

Define  wall  motion  relative  to  adjacent  cell  zone?  yes 

Apply  a  rotational  velocity  to  this  wall?  no 

velocity  Magnitude  (m/s)  0 

x-Component  of  wall  Translation  1 

Y-Component  of  wall  Translation  0 

Z-Component  of  wall  Translation  0 

Define  wall  velocity  components?  no 

x-component  of  wall  Translation  (m/s)  0 

Y-component  of  wall  Translation  (m/s)  0 

Z-Component  of  Wall  Translation  (m/s)  0 

internal  Emissivity  1 

External  Emissivity  1 

External  Radiation  Temperature  (k)  300 

Radiation  BC  Type  3 

X-Component  of  Radiation  Direction  1 

Y-Component  of  Radiation  Direction  0 

Z-Component  of  Radiation  Direction  0 

Theta  width  of  Beam  (deg) 

9.9923854e-07 

Phi  Width  of  Beam  (deg) 

9.9923854e-07 


(((constant  .  0)  (profile  ))) 

(((constant  .  0)  (profile  ))) 

(1) 

Apply  Direct  Irradiation  Parallel  to  the  Beam?  no 

Use  Beam  Direction  from  Solar  Load  Model  Settings  no 

use  Direct  and  Diffuse  irradiation  from  solar  Load  Model  settings  no 
Rotation  Speed  (rad/s)  0 
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x-Position  of  Rotation -Axis  origin  (m)  0 
Y-Position  of  Rotation-Axis  Origin  (m)  0 
Z-Position  of  Rotation-Axis  Origin  (m)  0 
X-Component  of  Rotation-Axis  Direction  0 
Y-component  of  Rotation-Axis  Direction  0 
Z-component  of  Rotation-Axis  Direction  1 
X-component  of  shear  stress  (pascal)  0 
Y-component  of  shear  stress  (pascal)  0 
Z-component  of  shear  stress  (pascal)  0 
Fslip  constant  0 
Eslip  constant  0 
Surface  tension  gradient  (n/m-k)  0 
Specularity  Coefficient  0 


Solver  settings 


Equations 

Equation  Solved 


Flow  yes 

Turbulence  yes 

Energy  yes 

Discrete  Ordinates  yes 

Numeri  cs 

Numeric  Enabled 

Absolute  velocity  Formulation  yes 
Relaxation 

variable  Relaxation  Factor 


pressure  0.3 

Densi ty  1 

Body  Forces  1 

Momentum  0.7 

Turbulent  Kinetic  Energy  0.8 

Turbulent  Dissipation  Rate  0.8 

Turbulent  viscosity  1 

Energy  1 

Discrete  ordinates  1 


Linear  solver 


vari  able 


Solver  Termination  Residual  Reduction 

Type  Criterion  Tolerance 
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Pressure 

X-Momentum 

Y-Momentum 

z-Momentum 

Turbulent  Kinetic  Energy 
Turbulent  Dissipation  Rate 
Energy 

Discrete  Ordinates 
Pressure-Velocity  Coupling 
Parameter  Value 


Type  SIMPLE 

Discretization  Scheme 
Variable 


V-Cycle 

Flexible 

0.1 

0.1 

0.7 

Fl  exible 

0.1 

0.7 

Flexi  ble 

0.1 

0.7 

Flexible 

0.1 

0.7 

Flexible 

0.1 

0.7 

Flexi  ble 

0.1 

0.7 

Flexible 

0.1 

0.7 

Scheme 


Pressure 

Momentum 

Turbulent  Kinetic  Energy 
Turbulent  Dissipation  Rate 
Energy 

Discrete  Ordinates 
Solution  Limits 


Standard 

Second  Order  Upwind 
Second  Order  upwind 
second  order  Upwind 
Second  Order  upwind 
First  Order  upwind 


Quantity 


Limit 


Minimum  Absolute  Pressure  1 

Maximum  Absolute  Pressure  5e+10 

Minimum  Temperature  1 

Maximum  Temperature  5000 

Minimum  Turb.  Kinetic  Energy  le-14 

Minimum  Turb.  Dissipation  Rate  le-20 

Maximum  Turb.  viscosity  Ratio  100000 
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Pipe  0.0002  m  thickness 

1 

2 

3 

4 

5 

6 

mass-weighted  average  enthalpy  (J/kg) 

innerjnlet 

-74019.76 

-74019.75 

-74019.58 

-74019.33 

-74019.73 

-74019.73 

inner_outlet 

-74012.18 

-74009.27 

-74009.52 

-74007.24 

-74002.15 

-73996.27 

difference 

7.58 

10.48 

10.06 

12.09 

17.58 

23.46 

J/kg-m  (quarter  pipe) 

37.9 

52.4 

50.3 

60.45 

87.9 

117.3 

J/kg-m  (full  pipe) 

151.6 

209.6 

201.2 

241.8 

351.6 

469.2 

Pipe  0.0004  m  thickness 

1 

2 

3 

4 

5 

6 

mass-weighted  average  enthalpy  (J/kg) 

innerjnlet 

-74019.76 

-74019.75 

-74019.602 

-74019.42 

-74019.74 

-74019.73 

inner_outlet 

-74013.13 

-74011.54 

-74011.203 

-74009.65 

-74006.59 

-74005.29 

difference 

6.63 

8.21 

8.399 

9.77 

13.15 

14.44 

J/kg-m  (quarter  pipe) 

33.15 

41.05 

41.995 

48.85 

65.75 

72.2 

J/kg-m  (full  pipe) 

132.6 

164.2 

167.98 

195.4 

263 

288.8 
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Pipe  0.0002  m  thickness 

1 

2 

3 

4 

5 

6 

total  heat  transfer  rate  (W) 

innerjnlet 

-28190.90608 

-28190.90614 

-28190.90621 

-28105.24270 

-28105.24278 

-28105.24349 

inner_outlet 

28188.02189 

28186.90937 

28187.00132 

28100.49146 

28098.55725 

28096.33527 

net  results 

2.88419 

3.99677 

3.90489 

4.75124 

6.68553 

8.90822 

W/m  (quarter  pipe) 

14.42095559 

19.98385586 

19.52443345 

23.75621690 

33.42766325 

44.54110230 

W/m  (full  pipe) 

57.68382234 

79.93542344 

78.09773381 

95.02486762 

133.71065301 

178.16440919 

Energy  (1  second,  J) 

57.68382234 

79.93542344 

78.09773381 

95.02486762 

133.71065301 

178.16440919 

mice  produced  (kg) 

0.00017271 

0.00023933 

0.00023383 

0.00028451 

0.00040033 

0.00053343 

AAbased  on  solidification  rate  of  ice 

volume 

1.88503E-07 

2.61218E-07 

2.55212E-07 

3.10528E-07 

4.36947E-07 

5.82216E-07 

AAbased  on  density 

Pipe  0.0004  m  thickness 

1 

2 

3 

4 

5 

6 

total  heat  transfer  rate  (W) 

innerjnlet 

-28190.906050 

-28246.585924 

-28246.585510 

-28105.242668 

-28105.242300 

-28105.242464 

inner_outlet 

28188.369141 

28243.442441 

28243.321652 

28101.407572 

28100.246475 

28099.730722 

net  results 

2.53691 

3.14348 

3.26386 

3.83510 

4.99583 

5.51174 

W/m  (quarter  pipe) 

12.68454113 

15.71741158 

16.31928808 

19.17548190 

24.97912600 

27.55871181 

W/m  (full  pipe) 

50.73816452 

62.86964633 

65.27715232 

76.70192759 

99.91650399 

110.23484724 

Energy  (1  second,  J) 

50.73816452 

62.86964633 

65.27715232 

76.70192759 

99.91650399 

110.23484724 

m,ce  produced  (kg) 

0.00015191 

0.00018823 

0.00019544 

0.00022965 

0.00029915 

0.00033004 

AAbased  on  solidification  rate  of  ice 

volume 

1.65805E-07 

2.05449E-07 

2.13316E-07 

2.50651E-07 

3.26513E-07 

3.60232E-07 

''''based  on  density 
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